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A quantum cloning machine is introduced that yields M identical optimal clones from N replicas 
of a coherent state and N' replicas of its phase conjugate. It also optimally produces M' = M + 
N' — N phase-conjugated clones at no cost. For well chosen input asymmetries N' /(N + N'), this 
machine is shown to provide better cloning fidelities than the standard (N + N') — > M doner. The 
special cases of the optimal balanced doner (N — N') and the optimal measurement [M — oo) are 
investigated. 
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The concept of cloning plays a central role in quantum 
information theory. It is for example crucial in quantum 
cryptography, as the optimal duplication of a quantum 
state directly determines the security of a cryptosystem 
jl],|]], or in quantum estimation theory, as cloning pro- 
vides a constructive way to achieve a generalised mea- 
surement (POVM) |J. The optimal N-to-M cloning 
transformations, which produce M clones from N orig- 
inals, have been found for both quantum bits 0] and 
continuous variables []5],[|. Interestingly, producing in- 
finitely many clones from N identical replicas of a quan- 
tum state is equivalent to performing an optimal mea- 
surement, which reflects the existence of a close link be- 
tween cloning and measurement theories. In the context 
of measurement, recent work has revealed that pairs of 
antiparallel qubits are intrinsically more informative than 
pairs of parallel qubits (j, a result that has later been 
extended to continuous variables: more information can 
be encoded in a pair of phase-conjugated coherent states 
\ip)\ip*) than in two identical replicas 1^)1^) @]< This 
property suggests that cloning machines with antiparal- 
lel input qubits (or phase-conjugated input modes) might 
yield better fidelities than standard N — > M cloning ma- 
chines, thereby opening a new avenue in the investigation 
of quantum cloning. 

In this Letter, we will focus on quantum information 
carried by continuous variables, and seek for a cloning 
transformation that, taking as inputs N replicas of a co- 
herent state \ip) and N' replicas of its complex conjugate 
\ip*), produces M optimal clones of \ip). The resulting 
concept of phase-conjugated inputs (PCI) cloning ma- 
chines will turn out to be closely connected to that of 
the amplification of light, just as for standard cloning 
P,pi]pjq|. As a matter of fact, it can be decomposed as a 
sequence of beam-splitters, a single non-linear medium, 
and another sequence of beam-splitters, which is con- 
sistent with the Bloch-Messiah reduction theorem | pd| . 
We will start by deriving the optimal canonical trans- 
formation that acts on two modes in a coherent state 
with respective mean values aip and pip* (where a, (3 are 
real while ip is a c- number) , and generates a mode whose 
mean value is yip, where 7 is real. Remarkably, this trans- 



formation will be shown to have a structure similar to 
that of a conventional phase-insensitive phase-preserving 
amplifier as defined in JT^|, where both the signal and 
idle modes are used as inputs. After having derived this 
transformation, we will apply it to the case of integer a 2 , 
P 2 , and y 2 , and see how it can be supplemented with 
beam-splitters to provide a PCI cloning machine for con- 
tinuous variables. This machine will be shown to produce 
M' = M + N' — N additional phase-conjugated clones (or 
anticlones) . The quality of the clones and anticlones will 
be discussed in the case of a balanced doner (N = N'), as 
well as for arbitrary input asymmetries N' / (N+N 1 ). The 
related question of the optimal measurement (M = 00) 
will also be treated. To our knowledge, the PCI doner is 
the first example of a new quantum information-theoretic 
process for continuous variables, for which no discrete- 
variable analogue has been found yet. 

Let {ai} and {b{\ (i = 1 ... 3) respectively denote the 
input and output modes annihilation operators of the 
cloning transformation. The indices i — 1,2 respectively 
refer to the input and phase-conjugated input modes, 
while i = 3 refers to an auxiliary mode. In full general- 
ity, we are seeking for a linear canonical transformation 

bi = Mijdj + Lijd'j (i, j = 1 . . . 3), (1) 

that meets the three following requirements (the sum 
over repeated indices being implicit). Firstly, starting 
with modes a\ and ai with mean values (a±) = aip and 
(0-2) = Pip* , we require = yip. We will only consider 
the case I7I > \a\, since, otherwise, the problem becomes 
trivial: one would just have to attenuate the input coher- 
ent state \aip) with an unbalanced beam-splitter, yield- 
ing a coherent state of amplitude yip. To simplify the 
problem, we may assume that /3 = 1, which amounts to 
susbtitute ip for flip. Then, we have: 

aM\\ + L12 = 7, 

Mx2 + aL n = 0. (2) 

Secondly, this transformation must obey the commuta- 
tion rules [bi,bk] = and [bi,b\] = 6ik (fi = 1), that 
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Thirdly, the noise of the output mode b\ of this transfor- 
mation should be minimum. 

Before sketching our calculation, let us note that a 
further simplification comes from the fact the the anni- 
hilation operators are defined up to an arbitrary phase, 
so that a transformation — * e Vi ai and b\ — ► e lv b\ al- 
lows us to take M\j and L\j as real and positive. Since 
we focus on phase-insensitive transformation, minimiz- 
ing the noise amounts to minimizing the sole quantity 
(A61) 2 = + 6 f iM - {h){b\) @. Thus, using the 

fact that (Aai) 2 = 1/2 for a mode a.; in a coherent state, 
we need to minimize 



(Ab 1 ) 2 = ±(M lj M lj +L lj L lj ), 



(4) 



under the constraints Eqs. (||) and (||). Rather than solv- 
ing this full problem, we use here a common trick in con- 
strained extremization problems that consists in solving 
a simpler problem with weaker constraints (bearing in 
mind that taking weaker constraints can only yield bet- 
ter solutions) and then checking that the solution of this 
simpler problem is one of the full problem. Specifically, 
we minimize (A&i) 2 taking into account the only condi- 
tion MijMij — LijLij — 1. Taking Eq. (|2|) into account 
and introducing a Lagrange multiplier A, we minimize 
the quantity + (7 - aM n ) 2 + (1 + o?)L\ x + M\ 3 + 
L\ 3 + A (Mfi - (7 - aMu) 2 - (1 - o?)L\ x + M 2 3 -L 2 13 -l), 
with respect to M\\, Ln, Mi 3 and L13. Some alge- 
bra shows that this problem admits only one solution 
A/13 = -^13 = Ln = M12 = 0, that is, the auxiliary mode 
is unnecessary. The optimal transformation has then the 
same structure as that of a phase-insensitive amplifier of 
gain G. Restoring (5, we get 



bi 

t>2 



VG~la\ + VGa 2 , 



with 



G = 



-cry - 



P 2 



(5) 



(6) 



It can easily be checked that, for (3 = (or a = 0), 
Eq. (^) reduces to a phase-insensitive phase-preserving 
(or phase-conjugating) amplifier as defined in Jl^] , and 
can be used to carry out the N — > M cloning (or phase- 
conjugating) transformation described in Q (or ||). 

Let us now turn to the special case where a 2 , 1 and 
7 2 are integers (which we will denote respectively as 
N,N',M). The transformation Eq. (§) can be used as 
the central element of a PCI cloning machine, which is 
covariant for translations and rotations in phase space 
(see Fig. H|). Indeed, the following procedure can be used 



to produce M optimal clones of a coherent state \ip) from 

(i) Concentrate the N replicas of \tp) stored in the N 
modes {q} (I = . . . iV — 1) into a single mode 01, this 
results in a coherent state of amplitude ^/N tp. This oper- 
ation can be performed with a network of beam-splitters 
achieving a iV-mode Discrete Fourier Transform (DFT) 
II , which yields the mode 
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and N — 1 vacuum modes. Similarly, concentrate the 
N' replicas of stored in the N' modes {di} (I — 
. . . N' — 1), into a single mode a 2 in a coherent state of 
amplitude VWip*, with the help of a iV'-mode DFT: 
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(ii) Process the modes a± and a 2 into a "phase- 
conjugated inputs" amplifier (PCIA), resulting in modes 
bi and 6 2 as defined in Eqs. (||) and (|6|). 

(iii) Distribute the output b\ into M clones {c[} (I = 
0...M-1) with a M-mode DFT: 
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where (k = 1 . . . M — 1) denote M— 1 vacuum modes. 
It is readily verified that this procedure yields M clones 
of IV)- 

Interestingly, the amplitude b 2 of the other output of 
the PCIA has a mean value y/M'tp*, with 



N — N' = M — M'. 



(10) 



Therfore, it can be used to produce M' phase-conjugated 
clones (or anticlones) of {d[} (I = . . . M' — 1), using 
a M'-mode DFT: 
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where {wk} (k — 1 . . . M' — 1) denote M 1 — 1 vacuum 
modes. Clearly, this procedure is optimal to produce M 
clones since its central element, the PCIA, is optimal, 
and the beam-splitters are passive elements. In addi- 
tion, the M' anticlones that are produced at no cost are 
also optimal. Indeed, our transformation produces M 
optimal clones and M > N, and is symmetric with re- 
spect to the interchange of labels 1 and 2. So, if our ini- 
tial problem was to produce M' optimal anticlones with 
M' > N', we would find the same solution. Noting that 
that M > N ^=> M' > N' [see Eq. @], it therefore is 
clear that our transformation yields both optimal clones 
and optimal anticlones. Furthermore, since the PCIA is 
linear and phase-insensitive, the resulting PCI doner is 
covariant with respect to translations and rotations of the 
state to be copied: all coherent states are copied equally 
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well, and the cloning-induced noise is the same for all 
quadrature components. 

Using Eqs. (||)-(||) and (|il"|), the noise of the clones 
and anticlones can be written as 
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where the gain can be reexpressed as a function of the 
number of inputs and outputs, 



Vg = 



VW - y/NM 
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As expected, the variance of the output clones exceeds 
1/2, implying that the clones are not exactly in the 
coherent state \ip). Instead, they suffer from a ther- 
mal noise with a mean number of photons given by 
(nth) = (G — 1)/A/. In other words, their P-function 
lOl is a Gaussian distribution 



i 



K(n t h) 



,-|C-0r/(«*h> 



(14) 



rather than a Dirac distribution P(£, £*) — <S^(£ — ip). 

Consider now the balanced case (N — N' , M — M'), 
for which simple analytical expressions of the noise vari- 
ances can be obtained. Taking the limit a — > (3 in 
Eq. (g) and replacing a 2 by N and 7 2 by M yields 
G = (M + N) 2 /4MN, so that the error variances of 
the clones and anticlones are 



(AcJ) 2 = (AdJ) 



(M - N) 2 
4M 2 N 



(15) 



Note that this balanced doner is optimal among all 
PCI doners in the sense that it minimizes (nth) f° r 
fixed N + N' and M + M'. It is convenient to char- 
acterize the quality of cloning in terms of the fidelity 

fm^j^ = (ip\p c \ip) /\(i/j\i/j)\ 2 where p c denotes the state 

of the clones. Using Eq. (H), we get 
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Let us now compare the production of M clones from N 
replicas and N antireplicas to the production of M clones 
from 27V identical replicas. The variance and fidelity of 
the clones fej obtained by standard cloning are given by 
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Of course, in the trivial case where M = 2N, standard 
cloning can be achieved perfectly, while the balanced PCI 
doner yields an additional variance (nth) = 1/(16/V). 
However, whenever M > 2N + 1, the (jy) — > M balanced 



doner always yields a lower variance (hence a higher fi- 
delity) than the 2N — > M cloning machine. The balanced 
PCI doner is also better for the anticlones: more anti- 
clones are produced at no cost, and they have a better 
fidelity. Indeed, a standard 2N — + M cloning machine 
produces M — 2N anticlones of fidelity 2N/2N +1, which 
actually is the fidelity of an optimal measurement of 2N 
replicas of \ip). In contrast, a PCI doner produces M 
anticlones with a higher fidelity, as given by Eq. ( |l6| ) . In 
particular, for a measurement (M — oo), we see from 
Eqs. (|l^) and ( |l7| ) that the additional noise induced by a 
PCI doner is 1/4ZV, that is, one half oi the noise induced 
by a standard 2N — > oo doner (\/2N). This reflects 
that more information is encoded in N pairs of phase- 
conjugated replicas of a coherent state than in 2N iden- 
tical replicas, a result which was proven for N = 1 in 
||. More generally, in the unbalanced case (TV ^ N'), 
it can be shown that the optimal measurement results 
in a noise that is equal to that obtained by measuring 
(y/N + yfW') 2 identical replicas of the input, in the ab- 
sence of phase-conjugated inputs. 

We have shown that the balanced PCI doner results 
in better cloning fidelities than a standard doner. More 
generally, we may ask the following question: If we want 
to produce M clones of a coherent state \tp) from a fixed 
total number n of input modes, N of which being in the 
coherent state \if>) and N' of which being in the phase- 
conjugated state \ip*) , what is the asymmetry a — N'/n 
that minimizes the error variances of the clones ? 

From Eq. (||), we see that for fixed values of the total 
number of inputs n and number of outputs M, the gain G 
(and thus the thermal noise of the clones (nth)) depends 
only on a, and varies as 
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In Fig. (|2]), we have plotted y (nth) as a function of a 
for n — 8 and different values of M > n. (Of course, 
only rational values of a are relevant here, but the whole 
curve has been plotted for simplicity.) In the trivial case 
where M = n = 8, the minimum additional variance is 
of course zero, and is obtained for a = 0. The cloning 
transformation is then just the identity. However, when 
M > n + 1 , using phase-conjugated input modes yields 
lower variances than standard cloning if a is correctly 
chosen (the lowest variance is attained for a ^ 0). Re- 
markably, the value of a achieving the minimum variance 
is not equal to 1/2 for finite M, that is the optimal in- 
put partition contains more replicas than antireplicas. In 
the limit of large M, however, the number of antireplicas 
achieving the lowest variances tends ton/2, and the curve 
G(a) tends to a symmetric curve around a = 1/2. This 
symmetry is not surprising, since M = oo corresponds 
to a measurement p],|T4| and we expect that measuring 
the value of ip from N replicas of and N' replicas of 
is equivalent to starting from TV' replicas of \tp) and 
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N replicas of \ip*). So, we have found that the optimal 
measurement is achieved with balanced inputs (N = N'). 
Finally, in the case where a = 1, the transformation con- 
sists in producing M clones of \ip) from n replicas of \tp*)- 
This is just phase-conjugation, for which we know that 
the best strategy is to perform a measurement ||. The 
additional variance is therefore given by l/n, which docs 
not depend on M. This explains why the curves converge 
all to the same point at a = 1. 

In summary, we have derived a continuous-variable 
cloning transformation using phase-conjugated inputs. 
This transformation has been shown to be decompos- 
able in a sequence of beam-splitters, a central ampli- 
fication stage, and another sequence of beam-splitters. 
A possible way to implement this central stage would 
be to use four-wave mixing. Two weak fields enter- 
ing the x' 3 ^ medium would then play the role of the 
phase-conjugated inputs, and energy would be brought 
to the system by two external modes in a large coher- 
ent state (see |13| for details). We have shown evidence 
that PCI cloning transformations outperform standard 
cloning transformations (taking only identical inputs) if 
the goal is to produce clones and anticlones of a state 
or to get knowledge about a state through measure- 
ment. The special case of the balanced doner, which 
produces M pairs of phase-conjugated clones from N 
pairs of phase-conjugated replicas, has been analyzed and 
was shown to be optimal. As far as we know, no qubit 
analogue of our transformation has been proposed yet, 
though it is very plausible that a cloning machine that 
produces M clones from N qubits in an arbitrary state 
and N' qubits in its orthogonal state can be defined. An- 
other possible extension of this work would be to study 
the case where the number of anticlones is a free parame- 
ter (in the PCI doner derived here, it is constrained by N, 
N', and M). Also, an interesting generalization would be 
to investigate fully asymmetric PCI cloning transforma- 
tions, that is, transformations whose output clones have 
all different fidelities. Finally, our work raises the ques- 
tion of why can a better cloning be achieved when phase- 
conjugated inputs are available, rather than identical in- 
puts. It is known that, in standard cloning, the spon- 
taneous emission occuring during the amplification stage 
is the physical mechanism that hinders perfect cloning. 
Thus, an open question left by our work is to understand 
why things happen as if comparatively less spontaneous 
emission occured in a PCI cloning machine. 
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FIG. 1. PCI doner that produces M clones and M' anti- 
clones from N replicas of l^} and N' replicas of \ip*). Modes 
are concentrated and distribitued by Discrete Fourier Trans- 
form (DFT). PCIA stands for a phase-conjugated inputs am- 
plifier. 
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FIG. 2. Cloning-induced noise standard deviation yj (nth) 
as a function of the fraction of the asymmetry a — N' /n, for 
n = 8 and several values of M/n. 
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